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Introduction 

The purpose of this interim report is to outline the methods being employed 
to meet the objectives of NASA Grant NAG-1-394 "Fundamental Studies of 
Structure Borne Noise for Advanced Turboprop Applications." The problem 
considered is the transmission of sound generated by wing mounted advanced 
turboprop engines into the cabin interior via structural paths and the 
comparison of the interior noise level due to the airborne and structure borne 
components. The goal is the assessment of the relative importance of the two 
sources of interior noise. 

The structural model employed is a beam representation of the wing box 
carried into the fuselage via a representative frame type of carry through 
structure. The structure for the cabin cavity is a stiffened shell of 
rectangular or cylindrical geometry. The structure is modelled using a finite 
element formulation and the acoustic cavity is modelled using an analytical 
representation appropriate for the geometry. The structural and acoustic 
models are coupled by the use of hard wall cavity modes for the interior and 
vacuum structural modes for the shell. The coupling method is similar to that 
employed by Craggs and Lau [1] and Dowell, Gorman, and Smith [2]. In the 
present study the coupling is accomplished using a combination of analytical 
and finite element models. The advantage is the substantial reduction in 
dimensionality achieved by modelling the interior analytically. 

Mathematical Model 

The mathematical model for the interior noise problem is readily 
demonstrated here with a simple plate/cavity system which has all of the 
features of the fuselage interior noise problem. The final model will replace 
the structural representation with a stiffened shell and wing. The final model 
for the cavity is no more detailed than described here, except that a 
cylindrical geometry can be used. 


A. Cavity Model 

The cavity is here modelled as rectangular with dimensions a, b, c. The 
face a,b is the plate, as shown in Figure 1. The acoustic environment within 
the cavity is governed by the field equations and boundary conditions 
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where p is the acoustic pressure, n is the outward normal at the cavity walls, 
and Cq are the ambient density and speed of sound in the air within the 
cavity, C is the portion of the cavity boundary which is rigid and S is the 
boundary of the cavity which is structural (the plate in the present case), w 
is the displacement of the structure at the fluid/structure interface, positive 
when directed outward. It is seen that the cavity acoustic pressure is driven 
by the acceleration of the structural walls. 

A weak formulation of the acoustic problem is written as 

ff/Wt Ptt3« -Jf * p O w tt^A 
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where the sum of the volume weighted residual of the field equation for a trial 
solution p from the class of continuous functions and the surface-weighted 
boundary residuals vanish for all weighting functions W-j from the class of 
continuous functions. If the weighting functions are constructed from the 
complete orthogonal set defined by solutions of 
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V 2 W * — W = 0 (5) 
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then the weak formulation can be written, after two applications of the 
divergence theorem, as 

Jff Wi [P t t + H 2 P]d¥ + ff p 0 c 0 2w i w tt dA = 0 
¥ S 

Next, the trial solution is constructed from the same complete orthogonal 
set defined by equation (5) 

P = E A i w i( x *y* z ) = C w C x ,y.2)] (A) (7) 

where []W(x,y, 2 )] is a row vector of acoustic hardwall mode shapes. Equation 
(6) now becomes 
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By virtue of the orthogonality properties of the acoustic hardwall mode shapes 
defined by equation (5) it is found that 

fffw T [ W ] dV- = C- N nn ^] (9) 

¥ 

where £ N nn <] is a diagonal matrix of "acoustic generalized masses. - Equation (8) 
becomes 


[' N nn .]{A tt } * Dn^C-N^fA) = -POC 0 2 Jf [W] T w tt dA do) 
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Equation (10) yields the unknown amplitude coefficients in the acoustic 
pressure representation of equation (7) for specified structural wall 
acceleration w tt . For the rectangular cavity discussed here the acoustic 
hardwall modes are 

W. = W. = cosHn — cosmir if- cosmir - (11) 

i Jlmn a d c 


B. Structural Model 

The structural model discussed here is a simple uniform plate of dimension 
a,b. The dynamics of the plate are described by the familiar result 

DV 4 W ♦ ph W tt = -(p a - p) (12) 


where 
D = 


Eh 3 
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E = Young's modulus 
h = plate thickness 
v = Poisson's Ratio 
p = plate mass density 


The applied pressure leading is p a and the pressure loading due to the cavity 
is p. The sign convention is set up so that the positive plate displacement is 
out of the cavity, consistent with the cavity model. Boundary conditions on 
the plate can be specified a; desired. For fuselage modelling the plate or 
shell model would have built in edges. For the present discussion, where it Is 



convenient to speak in terms of an analytical representation of the mode shapes, 
it is simpler to use a simply supported condition. 

A weighted residuals formulation for the plate can be written 

If *1 [°’ ,w ‘Ph w tt * Pa ' P3 dA = 0 03) 

S 

This becomes a classical Galerkin formulation if we choose the weighting 
functions ^ to be solutions of 

D - pho^ 2 ^ = 0 (14) 


with the same boundary conditions as the plate and also expand the solution w 
in terms of these functions 


w(x,y) =E^n( x »y) = l>Hq] (15) 

where [ ip] is a row matrix of eigenfunctions. 

If a finite element formulation is used the expansion is given by 

w(x,y) = [N]fr]{q} (16) 

where [N] is a global shape function matrix (interpolation matrix) and [y] is 
the model matrix (matrix with eigenvectors as columns) determined from a finite 
element solution of equation (14). The interpolation matrix is not explicitly 
defined globally, but consistent with finite element methods is defined on 
subdomains (elements). In the case of the plate the shape functions have 
continuous first derivatives. Equation (13) is written in terms of the unknown 
modal amplitudes as 
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The plate modes are orthogonal according to 
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Equation (17) now becomes 
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When finite element methods are used to obtain the mode shapes, equation 
(19) is written 

f' M nn')«lttJ * f<o 2 v]['M nn s](q} = ff [N) T pdA - [W ff [N] T P a dA ( 2 <» 
The coupled plate and cavity equations are now (10) and (20) 
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The integrals [v] T ff [N] T pd/ and PqCq 2 Jf [W] T w tt dA are C0U P lin 9 matrices 
which can be rewritten using equations (7) and (16) 

[Y] T fj [N] T pdA = [Y] T JJ [N] T [W]dA{ A) (22) 
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If we introduce the definition 
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then equations (21) can be written 
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where we have interpolated p a according to 

P a - (NllPa) 


with (p a ) being nodal values of p a (and derivatives if required) and 
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Equations (22) and (23) show that the coupling matrices are related to 
transposes of one another. The integrals in equations (22) and (23) and the 
integral defining [F] in equation (26) can be evaluated by usual finite element 
techniques of integration over subdomains (elements) using shape functions 
which are explicitly defined on the subdomain. 

In the k. ise of the plate, the shape functions are Hermitian polynomials 
which in the rectangular plate geometry load to conforming elements with 
complete slope continuity at element boundaries. In a completely consistent 
formulation of the integrals in equations (22), (23), and (26) the shape 
functions implied by the global shape matrix [N] should be Hermitian as 
required by the finite element analysis of equation (14) which leads to the 
modal expansion of equation (16). The finite element analysis of equation (14) 
to obtain the modal matrix [t] certainly requires slope continuity in the shape 
functions. The coupling matrix [0] does not require this much continuity and an 
interpolation based on only the nodal values of the weighting functions and the 
pressure or plate displacements is satisfactory. Hence, the most versatile 
procedure would evaluate the coupling matrix [6] and input generalized force 
distribution matrix [F] with the element shape matrices representing quadratic 
interpolation on 8 or 9 noded isoparametric elements. The modal matrix [y] is 
compressed to include only the displacement degrees of freedom. This procedure 
results in no detectable compromise in accuracy and allows experimentation with 
element types for the plate (or shell) model without requiring a reformulation 
of the coupling matrix computation. 

Equations (24) and (25) represent a set of coupled second order ordinary 
differential equations with generalized coordinates representing the plate 
modal amplitudes and the acoustic cavity modal amplitudes. Dissipation is not 
explicitly included in these equations but can easily be included as structural 
damping in the plate generalized stiffness matrix or as an equivalent viscous 
damping. Dissipative walls in the cavity can also be modelled. 

Equations (24) and (25) can be written in the form of a general dynamic 
system 

[M]{*x} * [C]{xJ ♦ [K]{x} = {Q} (27) 

where the generalized coordinates { x } are the vectors {q} and {A} and the 
matrices [M] and [K] can be constructed from equations (24) and (25). The 
viscous damping matrix can be constructed from equivalent damping for the 



structure or for dissipative walls. The generalized force vector { Q} is also 
easily identifiable in equations (24) and (25). 

The model for a stiffened fuselage structure is more complicated only in 
the structural modelling. All steps taken here will apply, however the finite 
element model for the structural mode shapes will be more complicated. 


Solution Methods 

It is presumed that the acoustic input to the structure is given in terms 
of its Fourier spectrum. That is 


Q(a») = J Q(t) e" ia)t dt 
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where Q(oj) is the Fourier spectrum in 
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Fourier Transform methods lead to the solution for the Fourier Spectra of the 
output generalized coordinates 


(X(gj)) = [D( 1 o)]' 1 {Q(u)} 


(30) 


where [D(oj)] is defined by 

[ D ( co) ] = [K] - w 2 [M] + ioi[C] (31) 

The acoustic pressure at any point within the cavity can be obtained from 
the acoustical modal amplitudes contained in { X ( co) } . Write 

{ A(co)} = [T a (co)]{Q(co)} (32) 

where [T a ( co) ] is the acoustic transfer function which is obtained by deleting 
rows from [D(co)]”* corresponding to the plate degrees of freedom. The acoustic 
pressure is obtained from 


P(x,y,z,&>) = [W(x,y,z)]{A(co)} 


= [W(x,y,z)]rT a (co)]{Q(a>)} 


(33) 



The energy density spectrum for the acoustic response at a point in the cavity 
’’s determined from the definition of the autocorrelation 

♦(T) - J P(t)p(t * T)dT 
• 00 

and its Fourier Transform 
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The inverse relationship 
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leads to the conclusion 
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This identifies the energy density spectrum |P(f)|^ as a measure of the system 
response since the "total energy" p 2 (t)dt is given by an integral over all 
frequencies of the energy density. From equation (33) 

|P(o .)| 2 = lQ*{w) } T [T a *(w)] T [W] T [W][T a ((i))]{Q(cD) J 

where the star represents the complex conjugate. The volume average of this is 

|P^l 2 * (Q*M) T [T a *M] T [ N nn ][T a (o.)]{Q(o.)J 

where [ N nn ] is defined by equation (9). Hence we write 

|PR ! 2 * (A*(u)} T [ N n „ ]{A(«)) 

*2> nn |A n(u) | 2 (34) 

A suitable measure of the cavity response is thus the weighted sum of the 
absolute values squared of the acoustic modal amplitudes. A suitable 
normalization of the cavity modes could be used to make the weighting factors 
N nn unity. 



Results 

The major effort to date has been expended in verification of the computer 
implementation. The following tasks have been undertaken: 

(1) Evaluation of plate elements with emphasis on convergence to exact 
free vibration mode shapes. This has established the Hermitian 
element as superior. 

(2) Evaluation of shell elements, again with emphasis on f-ee vibration 
convergence. The 32 degrees of freedom extension of the 
Bogner-Schmidt plate element produces the best results in the present 
case. 

(3) Detailed comparisons of fully analytic and combination finite 
element/analytic modes for the rectangular cavity - plate system. 

These comparisons were performed to verify the essential features of 
the coding relating to the coupling of the plate and cavity. Figure 2 
shows a response spectrum for a plate-cavity system generated from a 
finite element plate model a, id analytic cavity model. Comparison with 
computations of Craggs and Lau is good. Comparison with a fully 
analytic modal representation for both plate and cavity is excellent. 
The results for the fully analytic model is shown on a different scale 
in Figure 3. 

(4) Substantial progress has been made in the finite element modelling of 
stiffened plates and shells. These will be coupled with the analytic 
cavity models in the near future. 

Effort will be started to couple the wing to the plate and shell models and 
to synthesize a reasonable model for propeller interaction with the wing as a 
noise source model. An available propeller noise prediction program will be 
used to derive an equivalent airborne noise field for the purposes of comparing 
the airborne and structure borne interior noise. 
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Figure 2. Example of a Coupled Plate-Cavity Transfer Function. 
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Figure 3. Frequency response using analytic 
plate modes and cavity modes. 


